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ABSTRACT 


The RFSCAT(Rotating Fan-beam SCATterometer) is 
one of the two payloads of CFOSAT, the China-France 
Oceanography Satellite. It is a radar scatterometer designed 
to measure the electromagnetic back-scatter from wind 
roughened ocean surface. The operating frequency of the 
scatterometer is 13.256GHz (Ku-band) and has a swath 
about 1,000 kilometers. In this paper, based on the 
characteristics of echo signal of RFSCAT, the on-board 
processing of RFSCAT signal is introduced. A Doppler pre- 
compensation LUT and A slice division LUT are developed, 
and the signal processing algorithms are validated and the 
wind retrieval performances of RFSCAT are analyzed based 
on the well-developed data simulation system of CFOSAT 
RFSCAT. 

Index Terms— CFOSAT, RFSCAT, 
processing, simulation, wind retrieval. 


on-board 


1. INTRODUCTION 


The CFOSAT RFSCAT instrument uses two parallel, 
rotating slotted-waveguide array antennas with two fan 
beams corresponding to vertical and horizontal polarization 
that scanns in a conical pattern!'7!, The operating frequency 
of the scatterometer is 13.256GHz (Ku-band) and has a 
swath about 1,000 kilometer, which is centered on the 
spacecraft’s nadir subtrack and covering about 90% of 
Earth’s surface each day. The RFSCAT produces large 
overlaps on the observation regions within the total swath 
by successive sweeps and provides a large number of 
sigma0 acquisitions with diverse azimuth and elevation 
combinations for a single surface resolution cell. It is 
expected to acquire frequent global coverage and high 
quality wind acquisitions! 

Unlike QuikSCAT and ASCAT scatterometers, the 
Doppler effects of RFSCAT introduced by relative motion 
between the scatterometer and the earth with its rotation are 
3-D coupling distributions varying with latitude position of 
the satellite, as well as the azimuth and elevation of the 
resolution cells'!. Improper Doppler processing will cause 
serious received signal frequency deviations which will 
effect onboard data accumulation, slice division and the 
product quality of retrieved wind fields. 
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A frequency pre-shift LUT for transmitting signal by 
shifting central frequency of transmitting signal is used for 
narrowing the echo bandwidth and the tradeoff between the 
complexity of onboard signal processing implementation 
and accurate Doppler compensation is considered, and the 
residual Doppler frequency shift compensation and slice 
division processing onboard are combined into one step to 
implement via a well-designed 1-D LUT. The latitude and 
elevation dimensions are neglected and only the azimuth 
dimension is reserved by optimization in this processing. 

A data simulation system of CFOSAT RFSCAT is 
developed for signal processing algorithm validation and 
performance analysis. Unlike the end-to-end theoretical 
simulation model, which is independent of hardware system 
realization of scatterometer and suitable for the system 
design phase and theoretical performance analysis; the new 
simulation system considered the actual CROSAT RFSCAT 
onboard processing such as Doppler frequency shift 
compensation and slice division. Meanwhile, the orbit 
parameters and propagation are also well considered. The 
SGP4/SDP4 orbit propagation model is used for pulse and 
cell observation geometry computation. The input could be 
simulated random wind or ECMWF NWP wind. Taking the 
pseudo-L1B data as the main output simulation data, the 
CFOSAT RFSCAT data processing algorithms are 
developed and optimized. Based on the simulation system, 
the distance between ground slices after dimension 
reduction processing of slice are provided; and the Kpc 
coefficients of slices are computed and compared with 
approximate method. With the random input wind field, the 
1 revolution data of RFSCAT was simulated, and the wind 
speed bias and wind direction rms error are provided. 


2. ON-BOARD PROCESSING OF RFSCAT SIGNAL 


Doppler frequency effects of RFSCAT are induced by 
relative movement between scatterometer platform and 
earth and the rotation of the earth. Figure 1 is the Doppler 
frequency shift distribution and residual Doppler frequency 
shift distribution of CFOSAT RFSCAT, where the solid line 
denotes the beam far-end (outer part) and the dotted line 
denotes the beam near-end(Inner part) respectively, and two 
curves in the same color denote the same latitude interval. It 
can be shown that the Doppler frequency shift is significant 
determined jointly by azimuth angle, latitude and elevation 
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direction. In elevation direction, the maximum values in 
beam near end are about +200kHz, and the maximum values 
in beam far end are about +500kHz. In addition, for a 
specific azimuth angle, the Doppler shift difference between 
latitude intervals could be above 200kHz. 
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Figure 1: The Doppler frequency shift distribution and residual Doppler 
frequency shift distribution of CFOSAT RFSCAT (a) the Doppler frequency 
shift distribution of ascending orbit; (b) the Doppler frequency shift distribution 
of descending orbit; (c) the residual Doppler frequency shift distribution of 
ascending orbit; (d) the residual Doppler frequency shift distribution of 
descending orbit. The solid lines denote the beam far-end; the dotted lines 
denote the beam near-end. The red lines denote the [-85, -75] degree latitude 
interval; the green lines denote the [-5, 5] degree latitude interval; the blue lines 
denote the [75, 85] degree latitude interval. 


Pre-compensation by shifting central frequency of 
transmitting signal is used for compensating the mean 
Doppler with reference to antenna azimuth direction. The 
pre-compensation values are the averaged Doppler in 
elevation and latitude dimensions with opposite sign(64 
azimuth angle bins are adopted). It can be shown in figure 1 
that the residual Doppler shift after pre-compensation is 
reduced to about [-200,200]kHz. 

The transmitted pulse of RFSCAT is linear-frequency 
modulated chirp signal, the echoes are deramped by mixing 
with a chirped reference signal after digital downconversion 
to intermediate frequency (IF)°°. To extract the range 
information, an FFT is performed on the filtered signal and 
a periodogram is implemented by applying a magnitude 
squared operation. The effective bandwidth of deramped 
signal is about 500 kHz. The sampling rate is less than 
1.5MHz, a 4096-FFT will be adopted to extract the range 
information from the deramped signal, and the final bins are 
summed into 40 energy slices for downloading. In the 
onboard processing chain, a Doppler compensation module 
is needed to control the frequency characteristic of the 
reference signal, so the center frequency of deramped signal 
will accord with the band-pass filter”. 


In order to extract the range information accurately and to 
implement appropriate accumulation processing of 4096 
points in a pulse according to the approximate equal slice 
size on the ground, the residual Doppler frequency shift 
should be compensation carefully. However, the onboard 
processing complexity will be increased dramatically if the 
residual Doppler frequency shift is compensated completely. 
Figure 2 is the 3-D LUT dimension reduction and slice 
division diagram. After Doppler frequency shift pre- 
compensation, the echo frequency per pulse is computed. In 
order to optimize receiving, the sign of chirp rate is changed 
for forward and backward azimuth scan (the sign of chirp 
rate is negative for [0, 180) degree azimuth angle, the sign is 
positive for [180, 360] degree azimuth angle). The 
frequencies of equal ground slice size are obtained by spline 
interpolation according to the corresponding relationship 
between echo frequency and elevation angles and 3-D slice 
division LUT is obtained with respect to azimuth angle, 
latitude and elevation direction. An optimal LUT dimension 
reduction processing is adopted, and then the slice division 
is conducted on pulse echo power. In the LUT dimension 
reduction processing, for a specific azimuth angle bin, the 
accumulation point boundary of 1-D LUT is determined to 
ensure all the points in 3-D LUT are included and the 
accumulation slope of 1-D LUT is adjusted accordingly 
based on all the latitude accumulation slopes. 
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Figure 2: 3-D LUT dimension reduction and slice division diagram 


The 3-D LUT dimension reduction and 1-D LUT-based 
slice division processing bring about the significant 
improvement of onboard processing complexity. On the 
other hand, because the existence of residual Doppler shift 
and the use of simplified 1-D LUT, the slice size on the 
ground couldn’t be equal completely, and usually not all the 
downloaded 40 points per pulse stand for the approximate 
equal size slices. Generally, the effective downloaded point 
numbers corresponding to the slice numbers are between 25 
and 35. Considering the low SNR and possible echo energy 
leakage, some marginal slices should be discarded for 
sigma0 extraction. The alterable resolution slices couldn’t 
be available for CFOSAT RFSCAT. 
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3. SIMULATION AND PROCESSING 
OF RFSCAT DATA 


A data simulation system of CFOSAT RFSCAT is 
developed for data processing algorithm validation and 
performance analysis. Unlike the end-to-end theoretical 
simulation model, which is independent of hardware system 
realization of scatterometer and suitable for the system 
design phase and theoretical performance analysis; the new 
simulation system considered the actual CFOSAT RFSCAT 
onboard processing such as Doppler shift compensation and 
slice division. SGP4/SDP4 orbit propagation model is used 
for generation of the platform status vector and FRSCAT 
observation geometry. After the positions of ground WVCs 
are determined, the observation geometry of every pulse is 
computed. The pre-compensation LUT and slice division 
LUT are read in and the positions of slices, the X-factor and 
Kpc coefficients are computed. After the pseudo level 1B 
data is produced, the sigma0s of WVCs are generated and 
the wind retrieval and performance analysis are conducted. 
The wind retrieval and performance analysis module could 
accept both the simulated random wind field and external 
wind field such as ECWMF wind as input data. The MLE 
wind retrieval and median-filter-based ambiguity removal 
algorithm are adopted in the RFSCAT simulation system. 

The orbit parameters and system parameters of CFOSAT 
RFSCAT are listed in table 1. In this paper, only 1 
revolution of 25km resolution simulated scatterometer 
pseudo level 1B data is simulated and used for performance 
analysis. In the simulation, the value of normalized standard 
deviation of calibration error Kpr is set to 0.2dB. NSCAT-2 
GMF is used for simulation and wind retrieval, and the 
normalized standard deviation of the model error is 
computed as 


_ (0.0139? - 0.0381W +0.2753)-3 5 W<15 () 
pm 7 
0.01 W>15 
Where, Wis wind speed (unit is m/s), r is the resolution 
(unit is km). 


Kpc defined as the normalized standard deviation of the 
echo is used to indicate the magnitude of the measurement 
variations due to the instrument thermal noise and radar 
signal fading effects. Kpc coefficients are functions of the 
instrument antenna pattern, illumination geometry, 
modulation format, transmit pulse modulation, and receiver 
filter characteristics. 

Figure 3 is the slice sizes on the ground for different 
azimuth angle. It is showed that most of the slice sizes on 
the ground approximately equals to 10km. The slice with 
bigger size usually lies on the position far from beam center, 
especially for azimuth angle is 0 degree or 180 degree, and 
the slice size could be greater than 15 km. In Figure 4, the 
Kpc coefficients of RFSCAT are provided and the compare 
between the Kpc coefficients derived from the analytical 


Table 1 The orbit parameters and system parameters of CFOSAT 
CFOSAT Orbit Parameters CFOSAT system parameters 


Semi major (km) 6891.984 | Transmit power (W) 120 
AEs Carrier frequency 

Eccentricity 0.00123 (GHz) 13.256 

Inclination (deg.) 97.464 Pulse width (ms) 1.35 
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Figure 3: Slice sizes on the ground for different azimuth angle. The dash 
dotted line is minimum slice size on the ground w.r.t 40 downloaded points; 
the dotted line is maximum slice size on the ground; the solid line is the 
averaged slice size on the ground w.r.t specific slice. For the convenience 
of averaging, not all the slices are showed. (a) The slice sizes of 0 degree 
azimuth angle. (b) The slice sizes of 90 degree azimuth angle. (c) The slice 
sizes of 180 degree azimuth angle. (d) The slice sizes of 270 degree 
azimuth angle. 
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Figure 4: The Kpc coefficients for different azimuth. The red line denotes 
the coefficient A; the green line denotes the coefficient B; and the blue line 
denotes the coefficient C. The dotted lines stand for the coefficients 
obtained by approximate formula; the solid lines stand for the coefficients 
obtained by analytical formula. (a) The Kpc coefficients of 0 degree 
azimuth angle. (b) The Kpc coefficients of 90 degree azimuth angle. (c) The 
Kpc coefficients of 180 degree azimuth angle. (d) The Kpc coefficients of 
270 degree azimuth angle. 
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formula and those derived from the approximate formula. 
The approximate formula used for Kpc computation is 


= (2) 
B=- 6) 


C=—_ (4) 


Where, B, is measurement bandwidth of slice, T, is 


transmit pulse width, To is the receiving window width. 


Form figure 4, we find that the values between the Kpc 
coefficients derived from the analytical and that derived 
from approximate is very close. For CFOSAT RFSCAT, in 
most situations, the Kpc coefficients could be computed by 
approximate formula, which dramatically improves the 
computing efficiency. 
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Figure5: The wind retrieval performance estimation of CFOSAT RFSCAT. 
(a) is the averaged wind speed bias on different wind speed interval; (b) is 
the averaged wind direction rms error on different wind speed interval. 


Figure 5 is the wind retrieval performance estimation of 
CFOSAT RFSCAT. From figure 5(a), we find that the wind 
speed bias is greater at nadir and outer part of the swath 
than that of inner part of swath. Generally, the wind speed 
bias increases with the increasing of wind speed. The 
maximum averaged wind bias is about 2.5m/s which 
corresponding to the [15, 24] m/s wind speed interval. 
Consistent with wind speed bias, the wind direction retrieval 
performance is also deteriorated at nadir and outer part of 
the swath. On the contrary, the wind direction rms error 
decreases with the increasing of wind speed. The averaged 
wind direction rms error on wind speed interval [8, 15) m/s 
and [15, 20) m/s are almost below 20 degree. For wind 
speed interval [0, 4] m/s, the wind direction rms error on 
inner WVCs is close to 40 degree, and the wind direction 
rms error of nadir WVCs is about 60 degree. 


4. SUMMARY AND CONCLUSION 
The Doppler effects of CFOSAT RFSCAT are 3-D 


coupling distributions. A frequency pre-shift LUT is used, 
and the residual Doppler frequency shift compensation and 


slice division processing onboard are combined into one 
step to implement via a well-designed 1-D LUT. The 
latitude and elevation dimensions are neglected and only the 
azimuth dimension is reserved in this processing. A data 
simulation system of CFOSAT RFSCAT considering the 
actual CFOSAT RFSCAT onboard processing such as 
Doppler shift compensation and slice division is developed 
for signal processing algorithm validation and performance 
analysis. Based on the simulation system, the slice size on 
the ground after dimension reduction processing of slice are 
proposed; and the Kpc coefficients of slices are computed 
and compared. With the random input wind field, the 1 
revolution data of RFSCAT was simulated, and the wind 
speed bias and wind direction rms error are provided. 

Slice sizes on the ground shows that the on-board 
Doppler compensation and slice division processing 
perform well for CFOSAT RFSCAT, which obtain optimal 
tradeoff between the demand of subsequent data processing 
and the complexity of on-board signal processing. For 
CFOSAT RFSCAT, in most situations, in order to improve 
the processing efficiency, the approximate Kpc coefficients 
could be used, rather than analytical formula. Based on the 
well-developed simulation system, the echo signal and 
downloaded data of CFOSAT RFSCAT are obtained, and 
the wind retrieval performance is estimated. The statistics of 
wind speed bias and wind direction rms error show that the 
CFOSAT RFSCAT could provide high quality wind 
acquisitions. Furthermore, the wind retrieval performance 
could be improved by algorithm optimization in the future. 
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